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We construct all the unitary cubic curvature gravity theories built on the contractions 
of the Riemann tensor in D-dimensional (anti)-de Sitter spacetimes. Our construction is 
based on finding the equivalent quadratic action for the general cubic curvature theory and 
imposing ghost and tachyon freedom, which greatly simplifies the highly complicated problem 
of finding the propagator of cubic curvature theories in constant curvature backgrounds. To 
carry out the procedure we have also classified all the unitary quadratic models. We use our 
general results to study the recently found cubic curvature theories using different techniques 
and the string generated cubic curvature gravity model. We also study the scattering in 
critical gravity and give its cubic curvature extensions. 
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I. INTRODUCTION 



Recently, a new purely cubic curvature theory in (4 + 1) -dimensions with physically attractive 
properties was found in two separate works HI |2] . In constructing this theory, different guiding 
principles were used. For example, in [I] the new cubic theory was uniquely singled out from the 
requirement that, just like in the higher curvature part of the (2 + l)-dimensional new massive grav- 
ity (NMG) ^ , the trace of the field equations be second order and proportional to the Lagrangian 
itself, and that the field equations be second order for spherically symmetric solutions. On the 
other hand, the authors of [2] were searching for a cubic curvature theory with analytical, simple, 
spherically symmetric solutions and found the same theory, and coined it the "quasi-topological 
gravity". The theory becomes more interesting once lower powers of curvature are added to it. 
L)- dimensional extensions of the theory and its black hole solutions were found in [H El HI [5] . Also, 
extensive work on the model with regard to holography and the AdS/CFT appeared in [BHS]. Es- 
pecially in [71 [9], the trace anomaly structure of the dual conformal field theory is studied and a 
(holographic) c-theorem a la |10J and |11J was proposed for both even and odd dimensions. 

In this work, we will construct all the viable cubic curvature gravity theories with the guiding 
principle that the theory be unitarity around its (anti)-de Sitter [(A)dS] as well as its flat vacuum. 
We will find the most general cubic curvature theories that are perturbatively unitary in D dimen- 
sions. To check the unitarity of a given gravity theory around a fixed background, say, with metric 
gahi one has to find its propagator, that is the O [h?) expansion of the theory where h^h denotes 
the perturbation around the background. This is in general a highly complicated problem even 
for constant curvature backgrounds. But, the following observation somewhat simplifies the anal- 
ysis: In generic D dimensions, if a higher curvature gravity theory is to be unitary, its propagator 
should, necessarily, reduce to the propagator of one of the following four unitary gravity models 
(the details and the unitarity regions of these theories will be given in Section II). 

I. The cosmological Einstein theory, ^(i? — 2Ao), with a massless spin-2 degree of freedom 
which is unitary for any Aq and for k > 0. (We take the mostly plus sign convention.) 

II. Einstein plus Gauss-Bonnet (EGB) theory, \{R-2Aq) + 1 {Rlbcd - ^Kb + ^ ^it^i a 
massless spin-2 degree of freedom, but with an effective Newton's constant ^ = ^ + 

'^^D^){D^)^ where A is the effective cosmological constant whose reality puts a con- 
straint on the parameters. (Here and below, the effective Newton's constant refers to the 
coefficient in front of the linearized Einstein tensor.) For unitarity Kg > 0. 

III. Einstein plus scalar curvature theory, ^{R — 2Ao) -|- aR^, with a massless spin-2 degree of 
freedom and a massive spin-0 mode with an effective Newton's constant ^ = ^ + jy^Oi. 

For unitarity, > and the mass of the scalar mode = ^J^l^^^ (^£)^i^(j:)-2) should 

satisfy > in dS, and the Breitenlohner-Freedman (BE) bound |12 | 113] 

2 D-1 ^ 

> — 7^ tA, 

'-2(L>-2) ' 

in AdS. 

IV. Linear combination of the above three theories, ^{R — 2Ao) -|- aR^ + 7 {Rabcd ~ '^■^ab + -^^)' 
with a massless spin-2 degree of freedom and a massive spin-0 mode with an effective New- 
ton's constant = + jj^ct + ^^^^^^^^^^7- For unitarity, Kg > and the conditions on 
the mass of the scalar mode given in item III with the new Ke- 
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Note that we do not consider the purely quadratic theories, since we would like to have the well- 
tested Einstein's gravity in the infrared region. But, our analysis can simply be extended to such 
models by considering k — )■ oo. Of course, the above models are not unitary for generic values of 
the parameters (especially in a curved background) , there are constraints which we shall explore in 
more detail by looking at the one-particle exchange amplitude between two covariantly conserved 
sources [17]. We, also, do not consider the Fierz-Pauli massive gravity (with mass (hfiu ~ ^^))! 
even though it is unitary, it is not gauge invariant, so its propagator will not match any of the higher 
curvature models that we shall discuss. The classification of all the unitary cubic curvature theories 
that we shall present is exhaustive and devoid of any nonphysical or hard to justify assumptions. 
We have done a similar classification in three dimensions recently |I8J where we have found all the 
bulk and boundary cubic curvature gravity theories. 

The approach of analyzing the unitarity of the cubic curvature gravity theories by determining 
the quadratic curvature action which has the same propagator with the original cubic curvature 
action can also be used to construct the cubic curvature extensions of the critical gravity theory 
found in |14| [T5] . As in the case of the unitary theories listed above, one can analyze the propagator 
structure of the critical theory by taking the proper limits of the one-particle exchange amplitude 
given in [l7j and reveal the double pole structure of the critical gravity implying the logarithmic 
modes found in [231 - 125] . 

The layout of the paper is as follows: In Section II, we discuss the propagator and the spectra of 
quadratic curvature gravity theories in (A)dS and find the unitary regions that lead us to the above 
classification. In Section III, which is the bulk of this paper, we find the equivalent quadratic action 
that has the same free theory of a generic cubic curvature theory in D-dimensional (A)dS spacetime. 
In that section, we also discuss the unitarity of the previously found cubic curvature gravity theories 
mentioned above [H El HI IS] and the cubic curvature theory generated as an effective theory in the 
bosonic string theory [TH]. In Appendix A, we give the one-particle exchange amplitude of the 
critical gravity and construct its cubic curvature extensions. Some details of the computations 
regarding the equivalent quadratic actions are delegated to the Appendix B. 



II. PROPAGATOR STRUCTURE OF QUADRATIC CURVATURE GRAVITIES 

Any higher curvature theory built on the powers of the Riemann tensor and its contractions in 
the form / (Rabcd) (but not its derivatives) has a propagator that has the same structure as the 
propagator of the following most general quadratic theory 



I d'^x ^ [i (i? - 2Ao) + aR^ + fiRlb + 7 (Rlbcd " ^Rlb + 



(1) 



with possibly redefined parameters. Hence, to study the tree-level unitarity of a generic higher 
curvature theory, it suffices to understand the unitarity of ([T]): that is to find the constraints on 
the five parameters. In flat space with Aq = = A, it is well-known that for nonzero /?, there is 
a massive ghost and the theory is not unitary [2D] . (Note that Stelle's original result was only for 
L> = 4, in fact, it is now known that a nonzero (3 is allowed in D = 3 [3]. For D > 4, Stelle's result 
is still valid in flat space.) For curved backgrounds, the unitarity analysis is actually quite involved 
and was carried out in [T7], where the tree-level scattering amplitude between two background 
covariantly conserved sources was computed in the linearized version of the theory ([T]) augmented 
with a Fierz-Pauli (FP) mass term for the graviton. The FP mass helps to fix the gauge, and 
as long as the background is curved, after finding the propagator, it can be set to zero without 
introducing any van Dam-Veltman-Zakharov type discontinuity. 
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Since our subsequent discussions of the cubic curvature theories rest on the unitarity of the 
above quadratic theory, let us recapitulate the relevant analysis on the unitarity of ([T|. First of 
all, the two maximally symmetric vacua, with the Riemann tensor 



2A 

Rabcd — iy~(Z) 2y ^^'"^^^'^ dadQbcj i (2) 



of this theory satisfy 



One should keep in mind that 7 enters into the picture only for D > 5. The linearized field equations 
for the metric fluctuation hat = gab — dab around one of these constant curvature backgrounds 
become ED E2] 

Tab {h) = ^g^b + (2a + /3) (gabU - VaVb + J^dai^ R'^+P ipQab - J^dabR'-^ , (4) 

where we have defined the effective Newton's constant as 

1 1 4AD 4A „ 4A (L> - 3) (L> - 4) 

;^";: + z^" + z^^+ (^-i)(z)-2) ^- 

This expression was used in the list of the unitary theories given in the introduction. Here, Tab (^) 
contains the matter source Tab and all the higher order terms O . The linearized Einstein 

tensor, the Ricci tensor and the scalar curvature read 

Gab = ^ab - ^dabR^ " ^ _ ^Kb, (6) 



Rab = I {y^^ahbc + y^Vbhac " - V a^ bh) , i?^ = -Uh + V'^V'^Kb " ^^h. (7) 

It is important to realize that for all higher curvature gravity models, the building blocks in Q 
will remain the same. Only the coefficients will be affected. For example, the n*'^ order higher 
curvature terms such as r/i?", 77 (-^a?,)"''^) etc. will change ^ with an additional term proportional 
to T^A"^-*^, and will shift the coefficients (2a + j3) and /3 with a term proportional to r]A'^~'^. 

The tree-level scattering amplitude between two background covariantly conser ved, V„r"^ = 0, 
sources, as shown in Figure, is 



A = j d^x^gT'a,{x)h'^'{x). (8) 

Here, hab is the perturbation of the background metric created by one of the sources which is 
represented with Tab- Our normalization of the amplitude is in such a way that in four dimensions 
for K = SttGn we have the usual Einstein equation. The one-graviton scattering amplitude between 



Figure 1: Tree-level scattering amplitude between two background covariantly conserved sources via the 
exchange of a graviton 



the two covariantly conserved sources Tab and T', was calculated in [I7j^ as 
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(9) 



where c = i)^^2^^^ and A^^ is the Lichnerowicz operator acting on a symmetric rank-2 tensor. 
To get we have set the FP mass to zero as discussed above. To simplify the notation we have 
omitted the integral sign and the measure, and represented the Green's function as the inverse of 
the operator. Computation of (|9]) is a somewhat cumbersome problem, since one has to find the 
tensor Green's function corresponding to the equation O^j^'^'^hd = Tab where the operator C^^'^'^ 
has fourth order and second order derivatives, and a constant piece. The details are laid out in 
jl7| . therefore, we will not repeat it here, but just note that the best way seems to be to decompose 
hab into the transverse and traceless parts as 



Kb = hab + V(al4) + VaVbC/) + gab^fp, 



(10) 



then, gauge fix the equation with the help of a, say, FP mass term, and decompose Tab in a similar 
way. This allows one to identify the physical parts of hab, and relate them to their sources. As we 
noted above, unless D = 3, /3 introduces a ghost. This becomes more transparent if one organizes 
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the T — T' part of (|9|) and rewrites it as 



A 



--2T' 



ab 



{KePa + 1 
2 (/3 + c) 



A 



(2) 



4A 



+ 



c{D-l){D-2] 
8AD(3 
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r 
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-1 



rpab _|_ 



□ 



D-2 

-1 



-r 



Kepa + 1 □ + 



4A 
D-2 



-1 



(Ke/3n + l) (□ 



c{D-iy{D-2^ 
where the mass of the scalar excitation is 

1 2AD 
CK^ ~ {D-l){D-2 



T 

□ + 



1 



T 

(11) 



2AD 
{D -1){D- 2) 



n -1 



(12) 



If /3 7^ for generic D, the poles can be separated and one can "see" the ghost. Since we have 
already studied the D = 3 theories in detail [TH], we will consider D > 3. Therefore, we first set 
/? = 0, then the apparent pole □ = — (£)_^^(^_2) decouples. It is still a complicated task to explore 
the unitary regions of the remaining theory. One can somewhat simplify the analysis by relying 
on the unitarity of the pure cosmological Einstein-Hilbert theory whose unitarity has been checked 
with different techniques. One can write the one-particle exchange amplitude of the Einstein's 
gravity, that is a = 7 = 0, as 



2k 



Tab 



A 



(2) 



4A 
D-2 



rpob _j_ 



1 



D 



T' □ + 



4A 



D 



(13) 



which represents a unitary interaction via a massless spin-2 graviton for k > 0. This expression is 
our canonical example which will guide us understand the spectrum and the unitarity of the other 
theories. Note that this result reduces to that of ^23j in four dimensions. 

Among the unitary theories listed in the introduction, the EGB theory augmented with the 
term aB?, which was classified as the fourth theory, contains the first three cases in the proper 
limits, therefore, let us start with the scattering amplitude of that theory which reads 



2 k, 



rpl 



A 



(2) 



4A 
D-2 



rpab _j_ 



1 



{D-l){D-2 



1 

D-2 

1 



T' □ + 



4A 
D-2 



where the mass of the spin-0 mode given in (12) reduces to 

2AD 



D 



4{D- l)aK 
and the effective Newton's constant is 

1 



{D-l){D-2y 



1 AAD 4A iD -3)(D- 4) 

- H a H — -y. 

K D-2 {D-l){D-2) ^ 



(14) 



(15) 



(16) 



Let us now study in some detail the unitarity regions and exploit the constraints on the parameters 
of this theory for all values and signs of A. For A = 0, the amplitude reduces to the simple flat 
space limit 



A = 2k 



Kb {d'Y' T'^' + -^r [d^y'r 



{D-l){D-2 



-T' d 



D-2 



4 (£> - 1) Ka 



(17) 
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which is unitary for k > 0, and = ^^(^^Si^Ka ^ ^^'^ hence, a > 0. As expected in flat space, 
the Gauss-Bonnet term does not play any role at the linearized level, therefore, 7 is not restricted. 
On the other hand, in AdS (A < 0), the mass of the scalar excitation should satisfy the BF bound 

which allows negative mass squared values. There is a wide region of parameters that yield a 
unitary theory. To require attractive gravity or nonghost behavior, we impose Kg > 0, and we also 
demand that in the limit of vanishing A we have a unitary theory, and hence require «; > and 
a > 0. [Note that the latter two conditions can be removed to extend the unitarity regions.] All 
these conditions restrict 7 to 

\4kA D-2j {D-3){D-4)' ^ ^ 

In the dS (A > 0) space, there is no BF bound on the scalar mode, instead one has > 0, and 
the unitarity region reads 



7 > 



1 ^ aD{D-4:)\ {D-l){D-2) 
(D-2f J {D-3) {D-4)' 



4kA 



+ 



(20) 



Note that A is a function of Aq, a and 7. 

The scattering amplitude of the EGB theory, that is the second theory in our list, is obtained by 
setting a = in ( 14 ) and ( 15 ) which freezes the scalar mode and reduces the scattering amplitude 
to that of Einstein's gravity albeit with a modified gravitational constant. To get the i2— 2Ao-|-q;-R^ 
theory, one can simply set 7 



in (14) and (15) 



One can also find the scattering amplitude for the critical gravity by using ([9]) and the discussion 
on this is given in Appendix A. 



III. EQUIVALENT QUADRATIC ACTION 

In the previous section, we have classified all the unitary quadratic gravity theories in D- 
dimensional (A)dS spacetime. Our main task, which is the purpose of this work, is now to extend 
this analysis to cubic curvature theories. In D dimensions, the most general cubic curvature action 
constructed from the contractions of the Riemann tensor, but not its derivatives, reads [28^ 



xV-9 



^{R- 2Ao) + aR^ + (3Rlb + 7 (Rlbcd " ^Kb + 

K ^ 
+ (^iR^'^^^ Rp Rqust + a2R^'^^'^ Rpq^ Rrstu + Cl^R^'^ R^'^^pRrstq 

+ CL4^RRF'^^^ Rpqrs ~\~ CL^R^'^ R^^ Rprqs ~l" CLqRP'^ R!^)Rqr ~\~ Cl'jRRF'^Rpq + (jLgR 



(21) 



Altogether, we have 13 parameters and the main question is to find the ranges of these parameters 
for which the theory is perturbatively unitary around its (A)dS vacua (note that generically there 
will be 3 vacua). A direct way to answer this question is to compute the O (/i^f,) action and compare 
it with the unitary theories listed above. Computation of the O (/i^^) action directly by computing 
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the relevant tensors to O (Ji^i^ is a highly cumbersome analysis. To somewhat simplify the problem 
we will use a technique developed in [2^] and extensively used in [30j to study the unitarity of the 
Born-Infeld (BI) type gravities, and in [TBj to show the unitarity of the BI extension of NMG j31| 
and for the classification of all unitary cubic curvature gravities in three dimensions including the 
cubic curvature extension of NMG [32] • The technique consists of finding an equivalent quadratic 



action in curvature that has exactly the same propagator as (21). Moreover, as it will be clear 



below, this equivalent quadratic action also has the same vacua. This latter point is important since, 
with the help of the equivalent quadratic action, one avoids the derivation of the full nonlinear 
field equations to find the maximally symmetric vacua or determine the effective cosmological 



constant in terms of the 13 parameters of (21). The procedure of how one finds the equivalent 



quadratic action was described in p2j which we briefly repeat here. Consider a generic Lagrangian 
as a function of the Riemann tensor, C = y/—gF (^-R^d) • ^'"^ choose our independent field to 
be the up- up-down-down Riemann tensor which will simplify the computations, since this choice 
does not introduce the background metric or its inverse. In fact, in the Appendix B, to further 
simplify the computations, we will consider the scalar curvature R and the Ricci tensor to be 
the independent variables as well as the Riemann tensor R'^] Finding the vacua and the O (/i^) 
action means computing 



cd 



6C 



ab 



^gab 

gab 



^9 



ab 



(22) 



dab 

ab 



The first term in the right-hand side is irrelevant for our purposes, the second term gives the field 
equations and the third one gives the propagator. Finding an equivalent quadratic action means 
one finds a quadratic Lagrangian in curvature >Cquad-cquai = \/— 5 /quad-equal (^Rcd) 
same O {hP), O (h) and O {h?) expansions as the original Lagrangian C The equivalent quadratic 
Lagrangian can be obtained by expanding F [R!^ around the yet to be found constant curvature 
vacuum R^^ given in (^ up to second order as 



/quad-equal (^Rcdj = 



d'F 



L-^ {Red) J R'^b 



ab E>'*f' 

cd ^cd 



(23) 



Let us stress that \/— 5 /quad-equal (-?^cd) Same O {hP), O {h) and O (/i^) expansions as the 

original Lagrangian. This simple, yet remarkable, result follows from the fact that all the, say, 
O (/i^) terms of a given Lagrangian C = ^J—gF ^i?"^^ are in the form 



^o(ft,2) = 5g 



cd 



5C 



6g^'^5g 



ab 



59 



ab 



~9) 



O{h0) 



+ 



F 



O{h0) 



-9)om + 



0(hi) + Po(K^) 



~9) 



0(/l2) 



-9)o{h2)Fo{hO), 

(24) 
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where {^/—g)o{h") refers to the O (/i") expansion of \/—g, F terms follow similarly. Here, all the 
terms involved in 



Fo{hO) + Fo(h^) + Fo(h^) =F (^Rcd) + 

dF 



are given as 
dF 



cd 



0{h) 



+ 



dRfd 



r>ab 



+ 



_JPF_ 



^"5 

cd 



0(h) 



(25) 



As it is clear from the second line of this expression, the linear order in curvature also contributes 
to the O (h?) terms. Generically, all the terms up to order n in the expansion 



E 

i=l 



d'F 



^cd ~ ^cd 



(26) 



cd 



contribute to Fq^h")-^ 



Let us now apply these general ideas to (21). First of all, we do not need to consider the 



constant, the linear and the quadratic terms in the curvature, since they are already in the desired 
form. Let us rewrite the cubic action in such a way that R^^ appears as the only independent 
variable 



P {Red) =0'lRllR'^Rqu + «2-Rrf^pg-^iM + O^aRqRtpRll + a^RR^jR, 

+ a^R^R^gR'^^ + agRgR^R'^ + utRR^R'^ + a^R? . 



rs 
■pq 



"q "p "r 



"q^'-p 



(27) 



Here, note that we have chosen R^'^'^'^ R^ Rq^st specifically instead of R^'^'^^ R^ ^ Rqtsu which 
was given in [28j. Because it is not possible to put the latter in the up-up-down-down form, 
RPi'-^Rp^'^Rgtsu = EPlRgiRq\. But, the two terms are related to each other as 



1 



-It -iip ^ icqtsu — ^''p r ^''qust i . ^^pq ^'•rstu- 



(28) 



Leaving the details of the computation to the Appendix B, here we summarize the final results. 

(29) 



The equivalent quadratic action that has the same free theory as (21) becomes 



I 



1[R- 2Ao) + aR" + mlb + 7 (RlhrA - ^Rlh + R^ 



^ In fact, another way of finding tlie equivalent quadratic action is to observe that at the desired order of O (/i^ 
(Riem — Riem)'' = 0. We thank a referee for pointing this fact. 
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with parameters 



1 _1 

R K 

X 

k 



Ao 
a =a + 



12A2 
{D- lf{D-2f 

{D - 3) ai + 4a2 + 2 (D - 1) (ag + Da4) + (i? - 1)^ {a^ + ae + ^a? + ^^as 
DA ( ^ K 



K 3 (L/ — 2 
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2A 



{D-\)[D-2) 
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3ai - 602 - (-D^ - L> - 4j 04 + as + (-D - 1) (-03 + 2a7 + 'iDa^) 
2A 



7 =7 + 



{D -\){D- 2) 

X [-9ai + 24a2 + 4L»a3 + (2D - 3) 05 + (£> - 1) (4L>a4 + 3a6 + Da-j)\ 
2A 



[-3ai + 6a2 + (^ - 1) (as + Da,)\ . 



{D-l){D-2 

The equation that determines the effective cosmological constant is 
A-Ao 



2k 



+ /A^ = 0, 



/ = (Da + P 



{D-A) ^{D- 3) {D - 4) 
{D-2f '^^{D-l){D-2)' 



(30) 
(31) 

(32) 

(33) 
(34) 

(35) 



which, generically is a cubic order equation in A once k and / are explicitly written. To check the 
unitarity, one also needs the equivalent effective Newton's constant 



1 



1 



k^ D 



4AD 



-a + 



4A 
D - 1 



/3 + 



4A(L>-3) {D 



{D-l){D-2 



4), 
^7- 



(36) 



The equations from (30) to (36) can be used to classify all the unitary cubic curvature theories. By 



choosing the parameters of the equivalent quadratic theory, one can reduce it to any one of those 
four unitary theories listed in the introduction and studied in Sec. |llj For example, if a = ^ = 7 = 0, 
the free theory of the cubic curvature gravity reduces to the cosmological Einstein theory with a 
massless spin-2 degree of freedom. These three conditions (which are linearly independent when 
a, /3, and 7 do not vanish) reduce the 13 parameters of the cubic theory to 10 parameters. The 
unitarity condition k > gives a constraint on these 10 parameters yielding a wide range of unitary 
cubic curvature theories. Let us note that as long as any one of the lower order parameters ^, Aq, a, 
/3, and 7 do not vanish, the cubic equation (35) has at least one real root, meaning that the theory 
has a maximally symmetric vacuum with a nonzero curvature. Similarly, the cubic curvature 
theory can be reduced to the remaining four unitary theories, each time with only a couple of 
constraints on the parameters. This result classifies all the unitary cubic curvature theories in 
D dimensional (A)dS spacetimes. Therefore, it could also be used to check the unitarity of the 
existing cubic curvature gravity theories found with different guiding principles. Below we will 
study some examples such as the cubic curvature gravity obtained as an effective theory in the 
bosonic string theory [19], and the theories introduced in (TJ H] and [21 El E]- 



Using ( 30 36 ) , one can also find the cubic curvature extensions of the recently found critical 



gravity jl4t [TB] which is given in Appendix A. 
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A. The purely cubic theory 



In [H [2] , the following cubic curvature theory was introduced 

7 1 
6 



+ \R^miR,r - \RR^m„ + ^ii' ) ■ (37) 



The trace of the field equations coming from this action is proportional to the Lagrangian, and 
the action has just the first derivative of g"^^ (r) in the spherically symmetric ansatz. Note that 
there is no linear and quadratic terms in the curvature. One does not expect this homogeneous (in 
curvature) theory to have an (A)dS vacuum. This can easily be seen from the equivalent quadratic 
Lagrangian which follows from (29) 

ab\ _ A2 / 10A\ A / 2 , d2 



/quad-equal (-Red j ~ ~ ^ g" ) ^ \2 K^"''""^ ~ + ^ j " 



Using (35), one obtains A = which says that the equivalent quadratic action (38) vanishes. 
Therefore, there is no propagator in this theory and there are no local degrees of freedom, unless 
one introduces linear and quadratic or zeroth order terms in the curvature. (It is likely that these 
theories have local degrees of freedom at the nonlinear level as in the case of Chern-Simons theories 
which lack local degrees of freedom at the quadratic level, but have nonzero degrees of freedom 
at the nonlinear level |531.) In fact, for any purely cubic curvature theory, unless the trace of the 
field equations vanishes identically, one does not have local degrees of freedom. [Taking the risk of 
being too pedantic, let us note that this fact follows from the vanishing of O (h?) expansion of any 
cubic or higher order curvature invariant around flat space.] 

The generalization of (37) to D dimensions was given in [21 HI [U] in different forms. Here, we 
choose the action given in ^ 



Zo = a I dPx 



npqrs n t u „ , 1 



p r ^qtsu T _ 3) (^^ _ 4) 

3 (i^ - 2) RP'^R^%Rrstq + ^^^^'^^ RRP'^'^Rpqrs (39) 



3DRPm'''Rprqs + 6{D-2) RT'miRgr - ^^^^ ^^ RRPmpq + ^R^ 



In fact, a second action was also given in [2] 

RP^^^R^^'^Rrstu + (^2D-3)iD-4) ("^^ " 51? + s) R^^R^^Rrstq 
+ ^(^D^-4D + 2) RRP'^'^'Rpgrs + 12{D- 2) {D - 3) R^m^^Rprgs (40) 



3 
2 

+ 8 (Z) - 1) (D - 3) RPm^Rgr -6{D- 2f RRP^Rpq + ^[d^ -4D + 6j i?^ 

which is related with the previous one as Z'j^ = 2Zy) + where is the six-dimensional Euler 
density 

+ SRRP'^'^Rpqrs + 24RP'^R'''Rprqs + IGR^'^R^Rgr - URRPiRpg + R^. (41) 
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To find the equivalent quadratic action for we first use (28) to see the Lagrangian as a function 



of R'^. Then, from (|30 
theory 



31 



32 



33 



34), one arrives at the cosmological Einstein-Gauss-Bonnet 



/quad-equal (^cd) 



+ 



3{D-3) (3Z)2 - 15D + 16) 
2 {2D -3) [D-lf{D-2) 
3 (3L>2 _ 15D + 16) A 



a R 



2DA 



4{D -1){D -2) {2D -3) 



a 

abed 



3{D-2) 



4Rl + R^ 



(42) 



Again, one finds A = from (p5|), and the theory has no propagating degrees of freedom. 



The equivalent quadratic action for Z'j^ is proportional to (42) with a proportionality constant 

^^^3g2^i5D+i6 ^^'^ • Therefore, as expected in the purely cubic curvature theories, it, too, has no 
local degrees of freedom. On the other hand, if we add lower powers of curvature to or Z'^, 
breaking homogeneity, then we have theories with local degrees of freedom. As an example, let us 
consider Zd augmented with all the possible lower order terms 



-9 



1 



{R - 2Ao) + aR^ + /Ji?^^ + 7 (r^ 



.2 

abed 



^Rlb + R" 



+ Zd. 



(43) 



The equivalent quadratic action has the following parameters 
1 _ 1 3{D-3) {3D^ - 15D + 16) A^ 



2 (2D -3) {D- iy{D-2) 



Ao = -Ao + —j— -r 

K 3[D — 1) 



1 



a 



a. 



7 = 7 + 



3 (3L>2 _ 15D + 16) A 



(44) 



A{D- l){D-2) {2D -3) 

esult is not particularly illui 
a simple specific example. Consider Z) = 5, a = /3 = 7 = Ao = 0, then the three roots are given as 



where A is determined by (35), but the result is not particularly illuminating to depict. Let us give 



A = 0, A± = ±3\ 



(45) 



To have an (A)dS vacuum, na should be positive. For A-t, ^ = — | and Re = — f , therefore, to 
have Kg > 0, we must choose k < and a < 0. The theory is unitary for these values, but it has 
the wrong sign bare Newton's constant. Of course, around A = 0, the theory is unitary for k > 0, 
but the cubic terms do not contribute and the free theory is equivalent to the Einstein-Hilbert 
gravity. With nonvanishing Aq, the picture changes, there are regions where the theory is unitary 
for K > 0. Similar analysis can be done in D dimensions, but the punchline of this section is that 
the equivalent quadratic actions of Z£) and Z'j-^ are just the Einstein-Gauss-Bonnet theories. 



B. Cubic theory of Oliva Ray and Myers Sinha 



In [3] , a cubic curvature theory was introduced using the requirement that the trace of the field 
equations is a second order equation of the metric whose cubic curvature action in our notation 
reads 



I 



-{R- 2Ao) + 7 (Rlbed - 'iRlb + R^) + blW, + b2W2 + 63X6 



where W\;i are the inequivalent (for > 5) contractions of three Weyl tensors 

'a c ^bedf 



Wi = C'^'^C^'JCbedf W2 = C^^^'Cj^Ccdef, 



(46) 



(47) 
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and X6 was given in (41). The same theory was introduced using AdS/CFT and the existence of 



a simple holographic c-function in [9]. To find the equivalent quadratic action of this theory and 
study the unitarity, there are two ways. The first way is to express Wi and W2 in terms of the 
tensor structures that appear in (27) and use (29). We have carried out this lengthy computation, 



but since the result is the same as a simpler computation that we shall lay out now, we will not 
present it here. The second, relatively simple, way is to realize that Wi and W2 do not contribute 
to the propagator in (A)dS or any conformally flat space as a matter of fact, and therefore, from 
(|30l |3lJ [32] [33) |34]) one can find the equivalent quadratic action as 



/quad-equal ^-^cd) 



2Ao 8D{D-3){D-4){D-5)A^^ 
n {D 



lf{D 



+ 



+ 



1 12 {D- 3) {D 



4) {D-h)k\ 



7 + 



6(1) 



{D- 
4) {D 



2) 



R 



(48) 



5) A, 



{B-\){D-2) 



2 

abed 



Therefore, the free theory of (46) boils down to the EGB theory classified as the second unitary 



theory in the introduction. (The coefficient of the Gauss-Bonnet term can also vanish for a tuned 
value of 63 yielding an interesting theory that we deal below.) The effective cosmological constant 
satisfies 



A - Ao ^ {D - 3) {D 



4) 263A3 {D -3){D- 4) {D -5){D- 6) 



2k 



{D-l){D-2) 



+ 



{D 



lf{D 



0, 



(49) 



which has at least one real root for any value of the parameters in generic D > 6 dimensions. As 
expected for D < 6, Euler density xe does not contribute and for D < 4 the Gauss-Bonnet term 
does not contribute. The unitarity constraint 



1 _ 1 4A7 {D -3){D- 4) I263A2 (D -3){D- 4) (D - 5) {D 
Re {D-l){D-2) ^ {D-lf{D-2f 



6) 



>0, 



which somewhat simplifies to 

1 



3Ao 
kA 



2 2A7 (D - 3) (D - 4) 



> 0, 



(50) 



(51) 



c {D-l){D-2) 

restricts the parameters, but there is still a wide range of unitarity region. Let us specifically 



consider D = 7, and Aq = 7 = 0, for which one gets from (49) 



A = 0, 



A 



2^63 



(52) 



The existence of a constant curvature vacuum requires nb^ < 0. Assuming k > to have the usual 
Einstein-Hilbert gravity in IR, one should choose 63 < 0. However, this is in conflict with the UV 
unitarity, since Kg = — | < 0. Again, assuming nonzero Aq or 7 changes the picture and allows 



unitary regions which can be found from (49) and (51). Moreover, one can add aR and PR^f, to 
the action, and still have unitary regions. 



Let us now consider the case for which ( 46 ) reduces to the pure Einstein-Hilbert gravity at the 



free-level. For this case, one should tune the coefficient of xe (assuming D > 6) 

7(^-1) (^-2) 



6A(D-4) {D-5)' 



(53) 
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Then, (49) has two solutions 



3{D-l){D-2) 



l±Wl + 



16k;7Ao {D - 3)' 



3 (£» - 1) (L» - 2) 



(54) 



Further imposing the uniqueness of the vacuum A+ = A„, one has A = 2Ao and a tuned Gauss- 
Bonnet coefficient 



3(D-l)(L>-2) 
^ ~ 16kAo {D - 3f 

The theory has an equivalent effective Newton's constant 

4k {D- 3) 



> 0. 



(55) 



(56) 



Summing up the result: the following cubic action has a unique (A)dS vacuum around which it 
has a unitary massless spin-2 mode 



— (R — 2Ao) ^ ^-^ tJ- ( R^hrvi ~ 4-R^h + 



+ biWi + 62 W2 + 



{D-lf{D-2f 



64KAg {D - 3y {D -4){D- 5) 



X6 



(57) 



Furthermore, for Ag < 0, the theory is also unitary in the boundary of AdS, since the free theory 
is exactly like the cosmological Einstein's gravity with modified parameters: 



quad-equal 



D-6 
4k {D- 3) 



{R - 4Ao) 



(58) 



The fact that (57) and (58) have the same O {hP), O [h) and O (h?) expansions is quite remarkable, 
and is a testament to the effectiveness of the tools laid out in this paper in studying the spectra 
and the vacua of higher curvature gravity theories. 



C. Cubic gravity generated by string theory 

In [12], using the three and four point scattering amplitudes of strings, it was shown that 
heterotic and type-II superstring theories do not produce any R^ interactions, but the bosonic 
string has the following (in our basis) effective action 

J = i I d^x^ R+^ (rI,,, - 4Rl, + + (-2i?P^-i?;/i?,„,, + R^^^^R^I^Rrstu 

(59) 

where a' is the usual inverse string tension. In this action, the coefficients of i?^^ and R?' are 
not unambiguously determined by the computation of [TD], but chosen for their simplicity. In 
our analysis, we will keep this choice; but, without much effort, one can of course consider these 
coefficients to be arbitrary. Note also that there is no bare cosmological constant. Let us now 



consider (59) to be our microscopic theory, forgetting its string theory origin as an effective theory. 



and study its spectrum and unitarity. Clearly, around flat space, it describes a massless spin-2 
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unitary excitation. Around its assumed (A)dS vacua, the equivalent quadratic action follows from 



(|29l) 



quad-equal 



?,{D-lf{D-2f 



+ 1 + 



a'^K^ {D - 5) 



(D 



a' 



[R. 



2 

abed 



q'2A 



R 



+ 



2{D-l){D-2] 



[2R, 



2 

abed 



^ab 



■ (60) 



As before, A can be determined from the vacuum equation (35). Specifically, for the critical 
dimension of the bosonic string D = 26, we have A = 0, a'A ~ —2.37133, and a' A ~ 26025.2, but 
because of the last term in (60), the theory is not unitary in (A)dS. Even if one adds ^ 



k{D-1){D- 2) 



-g\R 



7 



-R 



ab 



(61) 



to the cubic action to turn the free part of the theory to the EGB theory with only a massless 
spin-2 excitation, the effective Newton's constant Kg becomes negative, meaning, one has a massless 
spin-2 ghost. Therefore, the theory is not unitary. 



IV. CONCLUSION 

We have systematically constructed all the perturbatively unitary cubic curvature gravity the- 
ories in D-dimensional (A)dS spacetimes. Our construction is based on finding an equivalent 
quadratic action that has the same vacua and the propagator as the original cubic theory. This 
equivalent quadratic action follows from a Taylor series expansion of the cubic action in powers of 



curvature up to second order. Using our general result, which was given in (29), we have studied 
the unitarity of the cubic curvature theories that appeared in the literature recently. These are 
the theories studied in [T| El SI El 13 E] that satisfy quite remarkable properties such as having 
second order equations and admitting simple c- functions. We have also studied the bosonic string 
generated cubic curvature gravity of [19 and found that it is nonunitary in (A)dS. We have also 
given unique vacuum cubic curvature extensions of the quadratic critical gravity theories. The tools 
we have presented here are very powerful in studying the free spectrum and the vacua of higher 
derivative gravity theories in (A)dS, and can be extended to any higher order. In fact, we have 
employed these tools in the infinite order Born-Infeld gravity theories |30J , and the Z)-dimensional 
Lovelock gravities [53] . 
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Appendix A: Scattering Amplitude and Cubic Curvature Extension of Critical Gravity 



Critical gravity |14| [T5j is defined by the action 



dPx^ 



-(i?-2Ao)+7CLd 

K 



(Al) 



A, and C, 



abed 



where the bare cosmological constant should be tuned to the value Aq = ^^^7(0- 3)^'' 
is the Weyl tensor. The spectrum of the theory involves only massless spin-2 degree of freedom. 
The effective Newton's constant is given as Kg = k ^1 — ^ j which has the opposite sign compared 
to the bare one. As we will see, the propagator of this theory has a double pole at the zero on-shell 
mass. 



1. Scattering in the critical gravity and its unitarity 

Scattering in critical gravity theory is a somewhat subtle issue, this is because critical theory 



was defined for T^b = 0, and the criticality strictly depends on Rl 



2A 
' D-2 



h = 0. If one wants to 



extend the critical theory to nonzero Tab, then, apparently, one should restrict it to the case when 
the trace of the energy-momentum tensor is zero, T = 0. For example, in four dimensions. Maxwell 
theory satisfies this condition and provides an example of an extension of the critical theory with 
sources. First of all, in the parametrization of ([T]), the critical theory corresponds to 



/3 



4a {D - 1) 
D ' 



A = An 



D 



a 



-fD {D - 3) 
{D-l){D-2)- 



Then, for T = 0, the scattering amplitude reads 



(A2) 



A = 2T' 



ab 



/?□ + 



1 



A 



(2) 



4A 



D 



(A3) 



(2) 

The Lichnerowicz operator AV acting on a traceless rank-2 tensor (pat yields 



-□ + 



4AD 



{D-l){D-2) 



and plugging the critical value of /3 in the form /3 



(D-l)(D-2) 

4AKe 



(A4) 



, one has the double pole structure 



A 



/3 



□ 



4A 



{D-l){D-2] 



□ 



4A 



{D-l){D-2] 



T 



•ab 



(A5) 



According to the old nomenclature, we have a degenerate Pais-Uhlenbeck oscillator; therefore, the 
two poles cannot be separated. With this observation, it is suggestive to think that the theory is 
free from the ghost, from which all the higher curvature theories (in > 3) with an ]r^^^^ term suffer, 
at the tree-level as long as /3 > (for AdS, this requires k < 0, and for dS k > 0). However, in flB], 
it is shown that logarithmic modes mix, in the Hilbert space, with the homogeneous (Einstein) 
modes; and if one tries to restrict the physical space to homogeneous modes only, the physical 
space will only involve vacuum state. 

Let us now consider the T 7^ case in the critical gravity. In this case, h cannot be zero, instead 
it is algebraically determined by the trace of the source 



h 



kT 

X' 



(A6) 
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which comes from Q. The degrees of freedom change and the scattering ampHtude is no longer 
given as ( A5 ) , instead one has 



A 



+ 



□ 



4A 



{D -1){D- 2) 



A 



(2) 



-r 



□ 



4A 



T' □ 



+ 



1 

A^(D-2) 
AD 

P{D-l){D-2) 



{D-l){D-2) 
4A 

{D-l){D-2) 



4A 
D-2 

□ + 



T 



ab 



4A 



D 



(A7) 



-1 



□ 



4A 

[D-l){D-2) 



□ + 



2kD 

{D-l){D-2) 



T. 



To see more explicitly the change in the degrees of freedom, let us compute the second order action 
in hab- The computation can be simplified by directly computing the linearized field equations 



ab 



4A 



iD-l){D-2] 



Gab + 



2%Ji) {aabU - VaVk - ^^9ab^ , (A8) 



and by integrating with the help of a reverse calculus of variations: — | / d^x ^J—gh"^^£ah where 
the overall factor comes from the self-adjointness of the involved operators. First, let us choose 
the gauge V^/iafe = Vf,/i as in IB], then = -^h. In this gauge, one has the linearized Ricci 
tensor 



^ab 



Va^bh - nhab + 



ADA 



{D -1){D- 2) 



h. 



4A 



ab 



and the linearized Einstein tensor 
1 
2 



-•L 
^ab 



^a^bh 



2A 



{D -1){D- 2) 
A{D-3) 



Qabh 



Uhab ) + _ 1) _ 2) + ^D-1){D- 2) 



(A9) 



(AlO) 



After using the gauge condition, one ends up with 

4A 



Tab 



□ 



{D-l){D-2) 



hab 



(All) 



/3A 



(D-l)(D-2) 



□ 



2A (D'^ -5D + 



Qabh. 



{D-l){D-2) 

After integration by parts, the second order action in this source-coupled critical theory becomes 

4A ^ 2 



'0(/l2) 



-9< 



2" r {D-l){D-2] 



h. 



ab 



2 ^{D-l){D-2) ^ {D-lf{D-2f 



(A12) 



which clearly shows that when h ^ 0, new degrees of freedom arise. To actually see these degrees 
of freedom and their masses, one has to do an orthogonal decomposition of the hab tensor as in 
(10), and further decompose h^^ into "spatial" tensor, "spatial" vector, and "scalar" parts in a 
suitable form of the (A)dS metric. This computation is not needed for the purpose of this work. 
(For D = 3, the computation was carried out for generic quadratic theory in [27].) 
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2. Cubic extensions of critical gravity 



Using the equivalent quadratic action for cubic curvature theories given in ( 29 34 ) , let us discuss 



the cubic curvature extensions of the critical theory found in [HI [T5]. To get a critical cubic 
curvature gravity in (A)dS, the parameters of the equivalent quadratic action should satisfy the 
following conditions. The first condition is 



Aa{D-l) 



D 



which removes the massive spin-0 mode, as long as 

K 



(A13) 



(A14) 



otherwise the trace of the linearized field equations is automatically zero, and therefore, Rl \s left 
undetermined. As a second condition, the mass of the spin-2 excitation is set to zero: 



+ 



4A/3 



{D -1){D- 2) 



0, 



which upon use of (A 13) yields 



1 . . ? 8Aa 
0=^ + ^^^+^- 



(A15) 



(A16) 



This second form makes the condition more transparent. In addition to these the vacuum equation 
has to be satisfied 



^ + /A^ 



f = [{D-l)iD-2)a + D{D-3)^] 



{D-A) 



D{D-1){D- 2) 



(A17) 



In the case of the quadratic critical gravity, for which f = f etc., (A16) defines the unique critical 



vacuum; therefore, one of the two roots of the second order vacuum equation is always noncritical. 
This could be remedied by imposing a unique critical vacuum which can be done in two ways: 
either one has coalescing roots or has / = 0. The first way contradicts (A16), but the second way 
leads to the unique vacuum critical theory 



4a (D 



D 



A = An 



D 



iKa 



a 



-fD{D-3) 
{D-l){D-2)- 



(A18) 



The cubic curvature extension is more subtle even though its equivalent parameters satisfy the 
same equations (A13), (A14), (A16), and (A17). The subtlety arises from the fact that tilde 



variables depend on A. Therefore, (A16) is a second order equation in A and apparently there are 



two critical vacua. Moreover, (A17) is a cubic equation in A and those two critical vacua should 



be roots of this equation. But, the third root will always be noncritical (assuming the first two 
are real). Therefore, just like in the case of the quadratic critical theory, one should impose the 
unique vacuum condition. This again can be done in two different ways: either one has coalescing 



roots or has a linear equation. Again the first way is in contradiction with (A16). The second way 
leads to A = Ao=Ao and k = k. As in the quadratic critical gravity, Aq is also tuned in terms of 
the other parameters. But, in this cubic critical gravity case, it satisfies a quadratic equation 



1 SApg _ 
D 



(A19) 
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where a is given in (32) (but, of course, two ofay's are eliminated in terms of the others). Therefore, 



for given parameters, say 7, ai, 02, 03, 04, 05, ag, there are two unique vacuum critical cubic 
curvature theories. For completeness, let us write the action 



1 



(i2-2Ao)+7CLd 



+ y aiR^'^^'^ Rp^j.^ Rqust + a2R^'^^^ Rpq^ Rrstu + O.-^R^'' R^'^^ pRrstq 

CL4^RRF'^^^ Rpqj-s ~\~ Oj^R^'^ R^^ Rprqs ~l" CLqR^'^ RpRqj- + d'jRRF'^Rpq + CLfi^R? 



(A20) 



with two of the parameters fixed as 
1 



a^ 



D{D-l){D-2) 



3 (£> - 6) ai - 24a2 - 4 (2L» - 3) 03 



AD {D - 1) a4 - {2D - 3) (L> - 2) 05 - 3 (D - 1) {D - 2) ae 



(A21) 



D^ {D-iy{D-2) 



2 [d"^ - SD + e) ai + 4 {fyD - 4) 02 + 2 (3D - 4) (D - 1) 



«3 



2L>2 {D - 1) 04 + (-D - 1) {D - 2Y as + 2 (D - 1)^ {D - 2) ag 



and Ao is a root of the following equation 

16kD {D - 3) [-3ai + 602 + {D-l) (03 + ^04)] 2 8/^7^ - 3) 



{D-lf{D-2f 



{D-l){D-2] 



(A22) 



Ao - D = 0, (A23) 



which gives a further restriction on the parameters coming from the reality of Aq. It is interesting 
to see that one can have 7 = 0, and generate a Einstein plus cubic curvature theory whose free 
theory is equivalent to the critical gravity, or one can have 7 7^ and extend the critical quadratic 
gravity to critical cubic gravity. 



Appendix B: Calculation of Equivalent Quadratic Lagrangian Densities 

To find the equivalent quadratic Lagrangian of a generic higher curvature gravity theory of the 
form C = ^—gF i^R'^ , one needs the following expansion 

/quad-equal (-^cd) (-^cd) + 



In 

- ^fr) (ci' - ^ri) • (Bi) 

ctfcm 
In 



dF 
dRt 



km 
In 



r>km 
^In 



d^F 



dR^!dR^ 
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Sometimes, computationally, it is easier to consider the Ricci tensor and the scalar as independent 
variables, and do the following expansion instead 



f ( J) T?"- T^o-b] T? ( T) Tja -nab 

/quad-equal y^i ' "-cd J ~^ I ^b ) ^cd 

OF 



+ 



dF 
dR 



R 



+ 



dF 



[r] - R] 

R 



+ 



+ 



d^F 
dR^ 



[r\ 

R 



km 75- 
In ~ ^ 



in J 



1 

+ 2 



J R 

d^F 



2 1 

R-R] +2 



d^F 

dR)dRf 



{r]-r]) [r^-r^ 

- R 



dRTzdR^^ 



( f}^^ _ nkm\ ( r>wy _ nwy 
y^^ln ^Hn J \^^xz ^^xz 



(B2) 



+ 



+ 



dF 
dRdm 



(i? - i?) (i?;- - R] 



dF 
IdRl 

dF 



dRdR\^ 



J R 



IX — It] \ rCi^ 



r>km\ 
^In ) 



dR^jdR'l^ 



y^j ~ ^j) \^ln ~ ^In ) ■ 



Of course (Bl) and (B2) give the same results, but we will use the latter below. Let us consider 



all the 8 cubic curvature terms in (21 ). 



1. Analysis of the three terms that do not involve the Riemann tensor 

Let us first find the equivalent quadratic Lagrangian for the simplest term F (R) = R^ which 
reads 



8D^A^ 12D^A^ 6DA , 

./quad-equal \it) — ,q . o it -f -r^ „it . 



(B3) 



{D-2f {D-2f" ' D-2' 

The term constructed from the Ricci tensor alone F = R^RpRf. has the following equivalent 
quadratic Lagrangian: 



Jquad-equalV-'T-ft/' ~ „^3 ' 7~i n ab- 



{D-2f {D-2f^^ ' D-2^^^'' ^-^^^ 

Finally, let us analyze the term involving the curvature scalar and the Ricci tensor, that is 
F (i?, i?^) = RR^Rp'. 



f . 8D'A' 12DA^ 4A , 2DA , 

/quad-equal [R, Rb) " _ " (^^^ + ^3^^ + ^^^-b- 



(B5) 



2. Separate analysis of the five terms involving Riemann tensor 



The term RP'^^R^p^ Rf^: 



Let F (i?^^) = RPIR'^'^R'^i. Then, the zeroth order part reads 

p { ^ab\ _ ^pq ^ru ^st _ SZJ {D - 3) 



(B6) 
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and the first order part reads 



dF 
dRt 



q TDPl T>rn 



12A2 



R 



{D - ly {D - 2Y 



{D - 2) 4C + 61,5^ 



(B7) 



In calculating the first order part, it is possible to take the derivatives with respect to the Riemann 
tensor such that the symmetries of the Riemann tensor are explicit. However, there is no need 
for such a care, because these symmetries will be taken into account upon multiplication with 

the diff'erence terms (^Rf^Y" ~ ^irT^ ^^'^ (-^f™ ~ -^fa") (^^xz ~ ^xf) which satisfy the required 
symmetries, and the parts that do not obey the symmetries give zero contribution. The second 
order part reads 



dF 



R 



3 (SyS^R^ + Rwk^y^m) 

( a' A^A^ A" — a' A^A^ A" 4- A^ a' A"A^ — A^A' A"A^ ^ 
\°y^k°w°m °y°k°m°w + °w°k°y'^m '^k^w'^y^m ) ■ 



{D-l){D-2) 

Therefore, the equivalent quadratic action for F (yR^ = Rf:jRp^Rq\^ becomes 

8D{D-3)A^ 12{D-3)A'^ 



/quad-equal (j^cdj 



{D-iy{D-2y {D-iy{D-2) 

6A „2 6A 



{D-l){D-2) 



iD-l)iD-2y 



(B8) 



(B9) 



b. The term RPlRf^RZ: 



read 



Let F = R^^R^pqRl^. Then, the zeroth order, the first order and the second order parts 



771 / T)ab\ -npq jytu r>rs ^ 

^ [Red ) - RriRpqRtu - 



{D-iy{D-2) 



3' 



(BIO) 



OF 
dRt 



24A2 



R 



{D-lf{D-2f 



(4C - ^A) , 



(Bll) 



dF 



dR^^dR^^ 



R 



{^k^m^wy + ^iv^y^kL) 



{D - 1) p - 2) (^^^-'^-^^" - ^kS^nSX + 5i5-5l5l - SlS^dl^si) , (B12) 



yielding an equivalent quadratic action 

32L>A3 



/quad-equal {Rccl) 



48A2 



R 



12A 



{D-iy{D-2y {D-iy{D-2y {d-i){d-2) 



-R: 



abed' 
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OF 



2m 



km /jj -|^\ /jj 2") yk"m "m'^k J ' 



R 



p-l)p-2) 



the equivalent quadratic Lagrangian becomes 



/quad-equal ^-^j ^cdj 



{D-l){D-2y {D-l){D-2y 



+ 



2DA 
D-2 



^abcd + 



8A 

{D-l){D-2) 



e. The term RlBLmi: 



Let F (-Rb,^^) = R^qRlR'^. Then, comput ing the relevant derivatives 



8DA^ 



{D-2) 



3 ' 



dF 



'JJR 



o a2 

{D-2f ' 



dF 
dRt 



A\2 



R 



(D - 2Y 



d^F 

dRpRf 



2R 



R 



{D-1) {D-2) 



dF 



dR^dR^:^ 



R 



the equivalent quadratic Lagrangian becomes 

/quad-equal {j^b i Rcdj ~ R 



2A 



2 2 {2D - 3) A 
(L>-2)=^ (D-2)2" ' p-l)p-2) ^ {D-l){D- 2) 
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